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The behavior of two-dimensional Ising spin glasses at the multicritical point on triangular and 
honeycomb lattices is investigated, with the help of finite-size scaling and conformal-invariance con- 
cepts. We use transfer-matrix methods on long strips to calculate domain-wall energies, uniform 
susceptibilities, and spin-spin correlation functions. Accurate estimates are provided for the location 
of the multicritical point on both lattices, which lend strong support to a conjecture recently ad- 
vanced by Takeda, Sasamoto, and Nishimori. Correlation functions are shown to obey rather strict 
conformal-invariance requirements, once suitable adaptations are made to account for geometric 
aspects of the transfer-matrix description of triangular and honeycomb lattices. The universality 
class of critical behavior upon crossing the ferro-paramagnetic phase boundary is probed, with the 
following estimates for the associated critical indices: v = 1.49(2), 7 = 2.71(4), r]i = 0.183(3), 
distinctly different from the percolation values. 

PACS numbers: 75.50.Lk, 05.50.+q 



I. INTRODUCTION 

In this paper we study two-dimensional Ising spin 
glasses, i.e., Ising spin-1/2 magnetic moments interact- 
ing via nearest-neighbor bonds of the same strength 
and random sign, drawn from a quenched probability dis- 
tribution: 

P(Jy) = p(5(Jy - Jo) + (1 - P) SiJ^l + Jo) . (1) 

Our interest focuses on the region of high (low) concen- 
tration p {1 — p) of ferro- (antiferro-)magnetic bonds, 
where even in two dimensions one can have order at 
T 7^ 0. A critical line on the T — p plane separates 
paramagnetic and ferromagnetic phases. Furthermore, 
for general space dimensionality d > 2 there is a second 
line of interest on the T — p plane, along which several 
exact results have been derived, known as the Nishimori 
line (NL)i*^. The shape of the NL is known exactly, and 
given by 
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A multicritical point is present, the Nishimori point 
(NP). The NP is believed^, to be located at the inter- 
section of the ferro-paramagnetic boundary with the NL. 
Many subsequent studies have taken this as a starting 
assumption, so far with consistent results, and we shall 
do so in the present work. As the shape of the phase 
boundary is known only approximately, e.g., from nu- 
merical studies, additional considerations are necessary 
if one intends to pinpoint the exact position of the NP. 

On a square lattice, a conjecture has been put for- 
ward^i^, to the effect that the NP should belong to a 
subspace of the T — p plane which is invariant under cer- 
tain duality transformations. For ± J Ising systems, the 
invariant subspace is given hy^: 

p log2(l + e-'-'o/^) + {l-p) log2(l + e^'"/^) = ^ • (3) 



Computing the intersection of Eqs. lO and (EJ, the exact 
location of the NP is predicted to be at p = 0.889972 • • • , 
T/Jo = 0.956729 ••• . This agrees well with earlier nu- 
merical estimates (though, in some cases, it is slightly 
outside estimated error bars). For detailed comparisons 
see, e.g., Ref. 0. 

^ Very recentljii, reasoning along the lines of Refs. 0, 
5 produced a conjectured duality relationship between 
locations of the NP on triangular and honeycomb lattices. 
By incorporating the NL condition, Eq. ((2l, considering 
lattices 1 and 2 dual of each other, invoking the replica 
method with n replicas and taking the quenched limit 
n ^ 0, and defining 



H{p) = -p logaP - (1 - p) log2(l - p) 



(4) 



it is predicted that, for mutually-dual systems with 
quenched randomness. 



H{pic) + H{p2c 



1 . 
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Using Monte Carlo simulations, the authors of Ref. 
established that Pc = 0.930(5) for the honeycomb, and 
0.835(5) for the triangular lattice. Using Eq. Q, these 
values imply that 0.981 < H (pic) + H {p2c) < 1.042, con- 
sistent with the conjecture Eq. ^ . 

Our goal here is twofold: first, to provide accurate 
checks of the location of the NP for both lattices, which 
will allow a more stringent test of Eq. lOl; and second, by 
invoking universality concepts, to gain more information 
on the universality class of the NP, through investigation 
of suitable critical properties on both lattices. 

Indeed, although many studies have dealt with the NP 
on square lattices, knowledge of the associated scaling in- 
dices is still restricted to (sometimes contradictory) nu- 
merical estimates. This is in contrast with the situation 
for pure discrete-symmetry systems in two dimensions, 
where it has been established that (i) all critical expo- 
nents are rational numbers belonging to a grid allowed 
by conformal invariance^, and (ii) for each universality 
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class the corresponding values have been unambiguously 
determined from the subset allowed by such grid, via ad- 
ditional exact results and/or numerical work. Even when 
(unfrustrated) disorder is introduced, significant progress 
can be achieved (for a recent review see, e.g., Ref. I^and 
references therein): though the connection to rational 
values of the exponents is lost, estimates obtained by 
various (analytical or simulational) methods are usually 
fairly consistent. 

Here we apply numerical transfer-matrix (TM) meth- 
ods to the spin-1/2 Ising spin glass, on strips of triangular 
(T) and honeycomb (HC) lattices of widths 4 < TV < 13 
sites (T) and 4 < iV < 16 sites (even values only, HC). 
In Sec.^ domain- wall energies are computed, and their 
finite-size scaling allows us to estimate both the location 
Pc of the NP along the NL, and the correlation-length 
index, yt = 1/v which governs the spread of ferromag- 
netic correlations upon crossing the ferro-paramagnetic 
phase boundary. In Sec. IIIIL uniform susceptibilites are 
calculated, and the associated exponent ratio 7/1^ is eval- 
uated. In Sec. IIVL we turn to probability distributions 
of spin-spin correlation functions, and their moments of 
assorted orders. These are shown to obey rather strict 
conformal-invariance requirements, once suitable adap- 
tations are made to account for geometric aspects of the 
TM description of T and HC lattices. Finally, in Sec.lvl 
concluding remarks are made. 
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Figure 1: Triangular lattice: domain- wall free energies, 
Eq. (O, along the NL (parametrized by p, see Eq. (O). Only 
data for odd are shown, in order to avoid cluttering. Error 
bars are of order of, or smaller, than symbol sizes. 



II. DOMAIN- WALL SCALING 

For pure two-dimensional systems, the duality be- 
tween correlation length ^ and interface tension a is well- 
establishedifi. For an infinite strip of width L, conformal 
invariance gives at criticalitjiii: 

LejL = TTT] , (6) 

where rj is the decay-of-correlations exponent, and ctl is 
the domain- wall free energy, i.e., the free energy per unit 
length, in units of T, of a seam along the full length of 
the strip: for Ising systems, gl^ ft- fli with /f (/f ) 
being the corresponding free energy for a strip with pe- 
riodic (antiperiodic) boundary conditions across. Finite- 
size scaling properties of ctl have been used in the study 
of critical properties of disordered systems as weliiSi, in- 
cluding an investigation of the NP on a square latticei^. 
With the above definition one has, for non-homogeneous 
coupHngs as is the case here, (jl = ^ln(A;^/Af) where 
In , In A^ are the largest Lyapunov exponents of the 
TM, respectively with periodic and antiperiodic bound- 
ary conditions across. 

We have calculated A^, Ag for strips of both T and 
HC lattices, usually of length M = 2 x 10^ columns, and 
widths N as listed in Sec. U (with the exception of iV = 16 
for HC). It must be recalled that both L in Eq. @ and 
the correlation length ^ (of which the surface tension is 
the dual) are actual physical distances, in lattice parame- 
ter unitsHiiLiSiil. Denoting by A the column-to-column 



distance by which the TM progresses in one iteration, the 
usual representations of the T lattice as a square lattice 
with a single diagonal bond, and of the HC as a "brick" 
lattice (i.e., with vertical bonds alternately missing), im- 
ply that L = C^, with Ct = 1, Chc = 3/2; At = \/3/2; 
Ai/c — \/3 (this latter is because two iterations of the 
TM are necessary in order to restore periodicity). The 
universal quantity rj is then given by 

= = / (2^/ V^^) (In - In kt) (T) 

TT |(7VV3/27r) (inA^ -lnAj4) (HC) . 

(7) 

For both lattices we scanned the NL, taking the respec- 
tive intervals quoted in Ref. Q as a starting guess for the 
location of the NP. 

For the T lattice, data for the scaled domain-wall en- 
ergy are shown in Fig. Standard finite-size scalingi^ 
suggests that the curves of Fig. would coincide when 
plotted against x = N^/" {p — pc). 

A quantitative measure of how good the data collapse 
is can be provided as follows. For trial values of {v,Pc) 
one calculates the per degree of freedom (x^d.o.f.) of 
a fit of the data to a phenomenological baseline curve 
f{x) (in the present case, since the curvature of data was 
monotonic, we found a parabolic form to be satisfactory). 
As the fractional uncertainties of data points were all of 
the same order, we used unweighted fits, i.e., the X^d o f 
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Figure 2: (Color online) Triangular lattice: scaling plot of 
domain- wall free energies, Eq. J^, against the finite-size scal- 
ing variable, N^^" (p — pc). The central estimates l/v = 0.67, 
Pc = 0.8355, have been used. Full line is quadratic fit to data, 
from which rj = 0.674(11) (see text). 



was calculated via 

X'd.o.f. = {Nd - M)-' J2 /(^«))' ' (8) 

where Nd stands for the number of data, A''^ — M is the 
number of degrees of freedom (M is the number of free 
parameters), {xi,yi) are the data points, and f{xi) are 
the values of the fitting function at the respective Xi . The 
use of uweighted fits is justified because all data used in 
each fit have similar fractional uncertainties. Therefore, 
the comparative analysis of different fitting parameters 
for a specified set of data will not suffer from distortions. 
This was the procedure used in all data collapse analyses 
in the present work. 

For domain- wall energies on the T lattice, we have 
found that the best collapse occurs for l/iy — 0.67(1), 
Pc = 0.8355(5). For the central estimates the X^doi 
is 3 X 10^^. Within the intervals of confidence given , 
the X^d.o.f. remains below 10~^. Fig. |2l illustrates the 
quality of plot obtained, when the central estimates just 
quoted are used. A parabolic fit to the scaled data gives 
r] = 0.674(11), where uncertainties in l/i/ and pc have 
been taken into account, in addition to those intrinsic to 
the fitting process for fixed values of these parameters. 

A similar line of analysis was followed for the HC lat- 
tice. Fig. El shows the unsealed domain- wall energy data, 
while Fig. 31 is a scaHng plot for the same data. The 
best collapse occurs for l/v = 0.67(1), Pc = 0.9325(5). 



Figure 3: Honeycomb lattice: domain-wall free energies, 
Eq. lO, along the NL (parametrized by p, see Eq. JSl). Error 
bars are of order of, or smaller, than symbol sizes. 

For the central estimates the x^d of is 7 x 10^^. Within 
the intervals of confidence given, the x^d.o.f. remains be- 
low 2 X 10~^. An estimate of r] from parabolic fits, 
with the same considerations used for the T lattice, gives 
77 = 0.678(15). 

The above estimates oipc for T and HC lattices, when 
plugged into Eq. Q, result in: 

77(pi,) + i7(p2c) = 1.002(3) . (9) 

This improves on the accuracy of the estimate given in 
Ref.Eby one order of magnitude, while still being com- 
patible with the prediction Eq. We view this agree- 
ment as a strong indication of plausibility of the conjec- 
ture exhibited in Ref. 0. 

As regards the correlation- length exponent, our esti- 
mate ly = 1.49(2) is incompatible with v = 1.33(3) quoted 
from the same sort of domain-wall scaling analysis ap- 
plied to the NP on a square lattice^^, but agrees well 
with ly — 1.50(3), found from mapping into a network 
model for disordered noninteracting fermions, via TMi^. 

Turning now to the exponent 77 given in Eq. ^ , it has 
been recalled, e.g., in Ref. 113, that in the presence of dis- 
order, the scaling indices of the disorder correlator (i.e., 
the interfacial tension) differ from those of its dual, the 
order correlator (namely, spin-spin correlations). Never- 
theless, the constraints of conformal invariance still hold, 
with the result that the amplitude of the domain wall 
energy remains a bona fide universal quantity^ft. For a 
square lattice, recent estimates give ry = 0.691(2) ^i?i^'?i^ft . 
This is sHghtly outside the error bars quoted here for the 
T lattice, but within the uncertainty given for HC data. 
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Figure 4: (Color online) Honeycomb lattice: scaling plot of 
domain- wall free energies, Eq. lO, against the finite-size scal- 
ing variable, N^^'^ (p — pc). The central estimates l/v = 0.67, 
Pc = 0.9325, have been used. Full line is quadratic fit to data, 
from which rj = 0.678(15) (see text). 



III. UNIFORM SUSCEPTIBILITIES 

We calculated uniform zero-field susceptibilities along 
the NL for both T and HC lattices, similarly to previous 
investigations on the square latticeSi. For the finite dif- 
ferences used in numerical differentiation, we used a field 
step 6h — 10""* in units of Jq. As in Sec^ we took the 
respective intervals quoted in Ref. Q as a starting guess 
for the location of the NP. 

Finite-size scaling argumentsi^ suggest a form 

X^ = Ar7/-/(Ari/-(p_pj) , (10) 

where xn is the finite-size susceptibility, and 7 is the 
susceptibility exponent. In order to reduce the number 
of fitting parameters, we kept l/iy and pc fixed at their 
central estimates obtained in Sec.^ and allowed j/f to 
vary. 

Within this framework, our best fit for the T lattice 
was for j/v = 1.795(20), as shown in Fig.^l For the cen- 
tral estimate the x^d of is 1 x lO"**. Within the intervals 
of confidence given, the x^d.o.f. remains below 3 x 10~^. 
These deviations are one and a half orders of magnitude 
larger than the corresponding ones for domain-wall scal- 
ing (see Sec. lTTji . 

We repeated the same steps for the HC lattice, with 
the results displayed in Fig. |£| The best fit now was for 
-f/iy — 1.80(4). For the central estimate the X^d.o.f. is 
1 X 10^^, two orders of magnitude larger than for the 





4 


_' ' 1 

: X 6 


1 1 1 1 1 


j/^ - 




3.5 






jT — 


X 




; A 10 


#11 J 




?^ 


3 


12 


• 


- 


1 












2.5 










2 


^tT I 


1 1 1 1 1 


, 1 ^ 



-0.02 0.02 
N'/^ (p-pj 

Figure 5: (Color online) Triangular lattice: scaling plot of 
uniform zero-field susceptibilities, Eq. I|in|l . The central esti- 
mates 1/v — 0.67, Pc — 0.8355, ^jv — 1.795, have been used. 
Full line is quadratic fit to data. 



T lattice. The lower quality of adjustment can be wit- 
nessed visually. Within the intervals of confidence given, 
the Q f remains below 3 x 10~^. Though the central 
estimates for the T and HC lattices are very close, the 
corresponding error bars differ by a factor of two. For 
the square lattice, we quote ^/v — 1.80(2)21, compatible 
with both values found here. 



IV. CORRELATION FUNCTIONS 

Our study of correlation functions is based on previous 
work for the square lattice^. We recall the following prop- 
erty, which has been shown to hold on the NL, for cor- 
relation functions dj between Ising spins en, cr j^i^i^^i^i? : 

(11) 

where angled brackets indicate the usual thermal average, 
square brackets stand for configurational averages over 
disorder, and I = 0,1,2,.... Denoting by P(Cij) the 
probability distribution function for the Cij, the pairing 
of successive odd and even moments predicted in Eq. lfTT|l 
implies that P'{Cij) = (1 — Cij) P{Cij) must be an even 
function of Cy, everywhere on the NL^. We have ex- 
plicitly checked that this constraint is obeyed by the dis- 
tributions generated for the T and HC lattices, within 
the same degree of accuracy as reported in Ref. for 
the square lattice. We shall not deal directly with the 
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Figure 6: (Color online) Honeycomb lattice: scaling plot of 
uniform zero-field susceptibilities, Eq. I|lf)ll . The central esti- 
mates 1/v = 0.67, Pc = 0.9325, ^/v = 1.80, have been used. 
Full line is quadratic fit to data. 



P{Cij) in what follows; instead, we concentrate on the 
scaling of their assorted moments m fe = [ Cj^] , especially 
in connection with their conformal-invariance properties. 

In contrast to the symmetry exhibited in Eq. lfTT|l 
which holds everywhere on the NL, conformal invariance 
is expected only where the NL crosses the phase bound- 
ary, i.e., at the NP. 

For pure Ising systems on a strip of width L of a square 
lattice, with periodic boundary conditions across, confor- 
mal invariance implies that at criticality, the correlation 
function between spins located respectively at the origin 
and at {x, y) behaves aA 



^ xy 



tt/L 



[sinh^(7ra;/L) + sin' 



1/2 



, = 1/4 . 
(12) 

For the T and HC lattices, the same is true, provided that 
the actual, i.e., geometric site coordinates along the strip 
are used in Eq. Ijl2|l . Thus, from the representation of the 
T lattice as a square (SQ) lattice with a single diagonal 
bond, and of the HC as a "brick" lattice, the respective 
SQ-like integer coordinates transform respectively 

into 



. .1 
— — z : V = 1 — - I 
2 2 



V3 



.7 + 1 



(T) ; 
(HC) , 



(13) 



where [X] denotes the largest integer contained in X. 



Figure 7: (Color online) Pure systems: ratio of calculated 
correlation functions at the critical point (c.f.) to asymp- 
totic form given in Eq. 11211 (cfO), against z = (sinh^(7ra::/L) + 
sm (Try/L))^''^ Main dia gram: T lattice (triangles), HC lat- 
tice (hexagons). Insert: SQ lattice. Strip width A'' = 10 sites, 
for all cases. 



Recall that L = N (T); L = iN/2 (HC), as ex- 
plained above. With R = {x^ + y^)^!'^ , the propor- 
tionality factor in Eq. Ijl2|l can be obtained from ex- 
act results (L,i? -> oo, i? < L), Cr = Ay^jR^I^, 
where = 0.703 38..., 0.668 65..., 0.767 04... re- 
spectively for X = SQ, T, HC24. Though strictly speak- 
ing Eq. II12|I is an asymptotic form, for the SQ lattice 
discrepancies are already very small at short distances^, 
and are even smaller for T and HC, as illustrated in 
Fig. H The horizontal axis in the Figure is the argu- 
ment z = (sinh2(7ra;/L)-Hsin^(7r2//L))i/2 of Eq. The 

range of z depicted corresponds to xjL < 0.6, i.e. (for 
strip width TV = 10 sites) up to, respectively, 5(HC), 
6(SQ), or 7 (T) full iterations of the TM. For larger x/L 
the angular dependence of z (through y) becomes less 
than one part in 10^. For z < 1 the discrepancy from 
Eq. lO is at most 0.4% for both T and HC, while in the 
worst case for SQ, namely {x,y)={l,l), it reaches 1.3%. 
For 1 < z < 3, the difference is < 0.1% for T and HC, 
and < 0.3% for SQ. 

The above analysis of conformal invariance of pure- 
system correlation functions indicates that, should sim- 
ilar trends hold at the NP of spin glasses, estimates of 
associated critical indices for T and HC lattices would be- 
have more smoothly than for SQ (since they rely on fits 
of numerically-calculated correlations to Eq. lfT2|l . with 
T] as an adjustable parameter). As in earlier work2>, we 
concentrate on short-distance correlations, i.e., where the 
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argument z is strongly influenced by y. Such a setup is es- 
pecially convenient in order to probe the angular depen- 
dence predicted in Eq. ltT2ll . which constitutes a rather 
stringent test of conformal invariance properties. 

We now turn to the quantitative analysis of the behav- 
ior of assorted moments mi of the correlation-function 
distribution, against z. Bearing in mind Eq. and 
following Ref. IQ, our goal is to extract the decay-of- 
correlations exponents r\2j+\^ via fits of our data to the 
form m2j+i ~ z~^^^+^. 

When one attempts such fits, several likely sources 
of uncertainty are present, on which we now comment. 
First, one has the finite width N of the strips used. In 
Ref.0, an extensive analysis of this point was undertaken, 
with the conclusion that, e.g. for TV — 10, finite-width 
effects are already essentially subsumed in the explicit 
L (i.e., N) dependence of Eq. lfT2|l . thus higher-order 
finite-size corrections most likely do not play a signifi- 
cant role. We shall assume that this is the case here as 
well, and restrict ourselves to = 10 for both T and 
HC lattices. Second, the finite length M of strips im- 
plies that averaged values will fiuctuate from sample to 
sample. Though the distribution itself (of, e.g., correla- 
tion functions) displays an intrinsic width which is a non- 
vanishing feature connected to the lack of self-averaging 
present at criticality, the average moments of the dis- 
tribution behave in the expected manner, namely, their 
sample-to-sample fiuctuations approach zero roughly as 
VM with increasing sample length MSSiSfe. Therefore, 
from a set of runs at assorted small values of M, one 
can infer what effect sample-to-sample fluctuations will 
have on results for larger M. In the calculation of re- 
sults shown below, we have used M — 10'', which implies 
a total of M' = 3.3 x 10® non-overlapping samples for 
our correlation-function statistics (because each sample 
needs three full iterations of the TM, in order to scan the 
set of lattice points of interest). For such value of M, the 
estimation procedure just outlined predicts fluctuations 
of order 0.1%, at most. 

Finally, one has the uncertainty in the location of the 
critical point. We have found that, in the present case, 
this is the main source of uncertainties for our data. 
Thus, e.g., with Pc — 0.8355(5) for the T lattice, averaged 
moments m2j+i taken at at the central estimate differ 
from those calculated at the edge of the error bar, by an 
amount increasing systematically with j, from ~ 0.7% 
for j = 0, to ~ 1.5% for j = 3. For HC, deviations follow 
the same trend against j but are slightly larger, ranging 
from < 1% for j = 0, to < 2% for j = 3. 

In Fig. IHlwe show data for the T lattice, taken at our 
central estimate for the location of the NP, p — 0.8355. 
The error bars, associated mainly to the uncertainty in 
Pc, as just discussed, are at most of order of the symbol 
sizes. Fig. El exhibits data for the HC lattice. Perti- 
nent comments are similar to those made above for the 
T lattice. 

In Table ffl we give numerical results of the fits illus- 
trated in Figs. IHl and El Though T and HC estimates are 
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Figure 8: (Color online) Triangular lattice: double- 
logarithmic plot of odd moments of the correlation- 
function distribution P(Cxy) against z = (sinh^(7ra::/L) -f 
sin^(7rj//L))^/^ (see Eqs. 11211 and ^^). Straight lines are un- 
weighted least-squares fits to data. Data taken at p = 0.8355 
for strip width A'^ — 10, M' = 3.3 x 10^ non-overlapping sam- 
ples in all cases. 



Table I: Estimates of exponents rj2j+i, from least-squares 
fits of averaged odd moments of correlation-function distri- 
butions. Data for A'^ = 10 and z < 1.6, assuming m2j+i ^ 
2-i2i+i rp. triangular lattice (this work); HC: Honeycomb 
lattice (this work); SQ: square lattice, calculated at the con- 
jectured exact location of the NP, see Eq. I^J (Ref. @). Last 
two columns: square lattice, authors as quoted. 



T 

0.181(1) 
0.251(1) 
0.297(2) 
0.330(2) 



HC 

0.181(1) 
0.252(1) 
0.296(2) 
0.329(3) 



SQ 

0.1854(17) 
0.2556(20) 
0.300(2) 
0.334(3) 



Ref. 

0.1854(19) 
0.2561(26) 
0.3015(30) 
0.3354(34) 



Ref. 

0.183(3) 
0.253(3) 



quite consistent with each other, and with the results of 
Ref. ITfll for j = and 1 both fall slightly below their SQ 
counterparts given in Refs. I^IT^ For j = 2 and 3, as a 
consequence of generally wider error bars, all estimates 
are broadly compatible with one another. 



V. DISCUSSION AND CONCLUSIONS 

We have used domain- wah scaling techniques in SecHTI 
to determine the location of the Nishimori point of Ising 
spin glasses on both the T and HC lattices. Probing 
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Figure 9: (Color online) Honeycomb lattice: double- 
logarithmic plot of odd moments of the correlation- 
function distribution P{Cxy) against z = {sinh? {-kx/ L) -f 
sin^(7rj//L))^/^ (see Eqs. 11211 and dJ). Straight lines are un- 
weighted least-squares fits to data. Data taken at p = 0.9325 
for strip width = 10, M' — 3.3 x lO'' non-overlapping sam- 
ples in all cases. 



Figure 10: Square lattice: semi-logarithmic plots of X^d.o.f 
for (unweighted) quadratic fits of domain-wall energies on 
SQ lattice to a dependence on the finite-size scaling variable 
N^^" [p — pc), against I /v. Each set of data corresponds to 
fixed Pc as fo llows: triangles, Pc = 0.889972 (conjectured ex- 
act, Refs. squares, Pc = 0.8894; hexagons, Pc = 0.8897; 
stars, jPe = 0.8903; crosses, Pc = 0.8906 (central estimate of 
Ref.Ea). Strip widths 4 < TV < 12. 



the temperature-concentration plane along the Nishi- 
mori line, we have obtained well-behaved curves of in- 
terfacial free energy; with the help of standard finite- 
size scaling techniques, we have extracted the estimates 
Pc — 0.8355(5) and Pc ~ 0.9325(5) respectively for the 
location of the Nishimori point on T and HC lattices. As 
a consequence of this, we have been able to refine the 
estimate of the quantity H{pic) + H{p2c) (see Eqs. Q 
and @), which has been conjectured in Ref. l^to be ex- 
actly unity. Indeed, our result given in Eq. Q is 1.002(3), 
which gives strong support to the conjecture cited. 

Furthermore, interfacial free energy data have allowed 
us to estimate the correlation-length exponent to be 
V = 1.49(2), in very good agreement with v = 1.50(3) 
from a mapping of the problem into a network model for 
disordered noninteracting fermionsi^, but incompatible 
with V = 1.33(3) from a TM treatment, presumably very 
similar to the present one, for the SQ latticei^. 

In order to investigate whether this latter disagreement 
might indicate a lattice-dependent breakdown of univer- 
sality, we calculated domain-wall free energies on the SQ 
lattice as well. Strip widths iV = 4 — 12 (both even and 
odd) were used, again with Af = 2 x 10® columns (ex- 
cept for iV = 12 where M = 1 x 10®). We scanned the 
region of the NL comprising 0.88 ^ p ^ 0.90, which in- 
cludes both the conjectured exact location of the NB^, 
namely Pc = 0.889972 • • • , and the estimate given in 



Ref.lil Pc = 0.8906(2). We found that scaled data col- 
lapse more smoothly for pc and v respectively close to 
0.889972 •• • and 1.5, rather than the values quoted in 
Ref. 0. This is illustrated in Fig. EH which exhibits the 
X^(j Q f for unweighted quadratic fits of scaled domain- 
wall energies to a dependence on the finite-size scaHng 
variable N^^" [p — pc)^ plotted against \/v. Each set 
of data corresponds to fixed Pc, see caption to the Fig- 
ure. We obtain u = 1.45(8), Pc = 0.8900(5), where the 
intervals of confidence given refiect the region in {v,Pc) 
parameter space in which the X^d.o.f remains below ~ 1.5 
times its overall minimum. Though the error bar for u 
is double that for T and HC lattices, the present esti- 
mate still encompasses the respective results for both, 
while excluding v = 1.33. For the domain-wall energy 
amplitude, we quote rj = 0.665(10), slightly lower than, 
but still compatible with, the values found in Sec. 
We conclude that our domain- wall energy data fully sup- 
port a picture of universal (i.e., lattice-independent) be- 
havior at the NP of T, HC, and SQ lattices. For all 
three lattices the correlation-length exponent is consis- 
tent with v = 1.50(3) of Ref.lTfll but most Hkely excludes 
ly = 1.33(3) of Ref.lil 

Our data for the uniform susceptibility, exhibited in 
Sec, mil do not scale as smoothly as the domain- wall en- 
ergies. Nevertheless, the application of finite-size scaling 
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ideas yields estimates for the exponent ratio ^jv which 
strongly support universal behavior at the NP, for T, 
HC, and SQ lattices. We recall that early work char- 
acterized the transition at the NP as compatible with 
the universality class of random percolation (see, e.g., 
Refs. HT^Ell for discussions of this point). In this con- 
text, we note that even our most accurate single result, 
namely ^/v = 1.795(20) for the T lattice, does not rule 
out the percolation valued {'-i/v)p = 43/24 = 1.7917 .... 
However, as explained below, consideration of the full set 
of results obtained here does support a scenario which 
rules out percolation-like behavior. 

Next, we turn to the investigation of correlation func- 
tions in Sec. IIVI The rapid convergence of T and HC 
results towards the asymptotic form, illustrated in Fig.d 
for pure systems, has translated to some extent into a 
discernible improvement on the accuracy of estimates for 
the disordered case. As explained above, for spin glasses 
on T and HC lattices the uncertainty in the location of 
the NP is the main source of fluctuations in numerically- 
calculated quantities. Thus, the relatively small uncer- 
tainties shown in Tablefflshow that the former effect com- 
pensates for the noise associated to the latter, at least 
partially. Compare, e.g., the T and HC columns with 
that for data taken at the conjectured exact location of 
the NP on SQi. 

The overall picture summarized in Table U clearly 
points towards universality of the several (multifrac- 
tal)^^'-'^ decay-of-correlation exponents. The small dis- 
crepancies observed, for j = and 1, between the T and 
HC estimates, and a subset of those obtained earlier for 
SQ, are likely to depend on details of the respective fltting 
procedures. One must note, however, that the j = and 
1 T and HC estimates are consistent with those derived 



m Ref.m This is similar to the case for the exponent 
V, in which our own result is compatible with the value 
found in Ref. 0, and not with that given in Ref. ^3 

Focusing now on j = 0, an unweighted average of all 
results of the corresponding line in Table ffl gives 771 = 
0.183(3). Considering the scaling relation ^/v — 2 ~ rji, 
one gets j/v = 1.817(3), which excludes {'j/v)p by a 
broad margin. Therefore, we quote the set of exponents 
1/ = 1.49(2), 7 = 2.71(4) rji = 0.183(3), distinctly differ- 
ent from the percolation values^l Up = 4/3, 7p = 43/18, 
Vp = 5/24. 

In summary, we have (i) produced accurate estimates 
of the location of the NP on T and HC lattices, which 
provide strong evidence in support of the conjecture ex- 
pressed in Eq. lOl; (ii) conflrmed that the critical proper- 
ties of the NP in two-dimensional systems are universal 
in the expected sense; and (iii) provided further evidence 
that such properties belong to a distinct universality class 
from that of percolation. 

As a final remark, we note that our discussion has been 
restricted to critical behavior upon crossing the ferro- 
paramagnetic phase boundary. The critical properties 
along the boundary Hue are of interest as welli2*21, and 
their investigation on the T and HC lattices would be a 
natural continuation of the present work. 
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